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Abstract 

Hadronic production of nonrelativistic boundstates of b' or isosinglet 
quarks D with suppressed weak decays is investigated for LHC and 
SSC energies. QCD corrections to production and decay rates are 
incorporated. Large rates for final states from r]y — > HZ are predicted. 
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The theoretical and experimental aspects of hadronic quarkonium production 
have attracted a great deal of interest after the original discovery of J/\I/ and 
T. Gluon fusion seemingly dominates the cross section with a small admixture 
of quark- ant iquark annihilation in pp and ir~p collisions at lower energies. Given 
sufficiently high center of mass energy and luminosity, the same parton processes 
could also lead to toponium production at a future hadron collider. The produc- 
tion rate of 1 bound states through the reaction g + g — > 1 + g is tiny. 
The cross section for gluon-gluon fusion into the 1 Sq state r] t , however, would 
be sufficiently large to produce a large number of resonances. The experimental 
signal that might conceivably be accessible under realistic assumptions is the 77 
decay mode of r\ t |1J. However, this mode is severely reduced by the tiny branch- 
ing ratio of order 10~ 4 down to 10 -5 — a consequence of the large single quark 
decay rate of 7] t . Detailed investigations which include higher order corrections 
for the production rate and a realistic QCD potential show that the window for 
the discovery of rj t is extremely narrow. Fairly optimistic assumptions have to be 
employed to access quark masses even in the range of 100 to 110 GeV 0, [|. 

However, the situation changes drastically if the single quark decay mode of 
the 1 5'o bound state is absent or at least strongly suppressed, as is the case for a b' 
with suppressed couplings to the lighter u, c, and t quarks, or for isosinglet quarks 
denoted in the following by D. The branching ratio for decay modes of these 
bound states with truly characteristic final states such as 77, jZ, ZZ, WW, HZ 
are no longer drastically suppressed and may even become dominant, as is the 
case for the HZ mode in the high mass region. These would be easily accessible in 
future experiments at LHC or SSC, for quark masses even up to 1000 GeV. Pre- 
dictions for these reactions in Born approximation have been derived in references 
1, 

In a recent paper 0, |), QCD corrections were evaluated for the production 
of nonrelativistic 1 S'o bound states in hadronic collisions. These corrections are 
positive and fairly large, justifying the renewed investigation of exotic bound state 
production. The corrections can be also applied to the production of bound states 
composed of a b' with J 3 = —1/2 and cq = —1/3 or of a D quark with J 3 — 
and cq = —1/3. 

In the following we specify the general structure of the cross section for r\yp 
production in hadronic collisions 

h 1 (P 1 ) + h 2 (P 2 )^r ]b , >D + X (1) 

in the framework of perturbative QCD. In (|I|) hi and h 2 denotes unpolarized 
hadrons with momenta P{. The hadronic cross section in NLO is thus given by 

a H (S) = J dx 1 dx 2 f^(x 1 ,Q 2 F )f^(x 2 ,Q 2 F )& ab {s = x 1 x 2 S,a s (fJ, 2 ), t J, 2 ,Q%) (2) 
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where one sums over a,b = q,q,g. fa(x,Q 2 F ) is the probability density to find 
parton a with fraction x in hadron h if it is probed at scale Q 2 F and a ab denotes 
the parton cross section for the process 



a(xiPi) + b(x 2 P2) -> Vb',D + X 



(3) 



from which collinear initial state singularities have been factorized out at a scale 
Q 2 F and implicitly included in the scale-dependent parton densities f^(x,Q F ). 
The partonic cross sections a ab in NLO for the gg, qg and qq initiated reactions 
are given in ||. They are proportional to the square of the bound state wave 
function at the origin R(0) and hence dependent on the QCD potential. For 
the evaluation we employ the two- loop QCD potential Vj ]/]]. The numerical 
results depend critically on the value of the QCD parameter A^,. The results 
for the dimensionless quantity \R(0)\ 2 / (M v = 2m b /,o) are displayed in fig. 1 

(200, 300 and 500 MeV) to indicate the characteristic 



for different values of A 
uncertainty in the predi 
and 0.127 respectively (and to A 
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uncertainty in the prediction. These A values correspond to & jjg = 0.109, 0.1165 



• (5) - 0.132, 0.207 and 0.366 MeV). For the 
parton distribution functions we employ MT set Bl with A^- = 194 MeV f| 
and work in the DIS factorization scheme. The renormalization and factorization 



scales are fixed at M = M v . We use the MS definition of a s at two-loop accuracy 
with six flavours 
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6(153- 19ra/) In ln(Q 2 /A 2 
(33 - 2n/) 2 ln(Q 2 /A 2 ) 



(4) 



and A^r = 58 MeV to be consistent with the A^- value from the parton distri- 
butions. 

The if— factor and the production cross section are shown in figs. 2 and 3a 
respectively as functions of the bound state mass for three different A values in the 
potential. The uncertainty from the choice of the parton distribution functions is 
shown in fig. 3b. The additional uncertainty from the choice of the renormalization 
and factorization scales is estimated to be less than 10% ||. 

To obtain predictions for the interesting reactions these production cross sec- 
tions have to be multiplied by the various possible branching ratios. It is con- 
venient to normalize the decay rates relative to the two photon decay rate T 77 
with 

Mo)| 2 



77 



12eja 2 



M 2 



(5) 



With R ab := r afe /r 77 (ab = -yZ, HZ, ZZ, WW, gg, tt) one finds 
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where 



1 1 /2p\ 3 

* w = 2e 4 Q (4sin 2 M 2 (l-p) 2 ImJ (9) 

^ = 9^ (10) 

R* = ^M^(^) (ID 



(12) 



Ml KM J 

Pl = M?/M 2 , p = 2 (1/4 - p Q , + Pw ) 



w Q = 2/ 3Q - 4e Q sin 2 % , a Q = 2/ 3 q y = 2 sin 26> w 

QCD corrections to the decay rates will be discussed below. In eq. fll2|) /3Q 
denotes the third component of the weak isospin of the Q— type quark and 6w is 
the Weinberg angle, p denotes the three-momentum of the final state particles. 
The single quark decay mode which dominates for rjt proceeds through the decay 
of the heavy quark into b+W. For a heavy sequential quark b' of isospin J 3 = —1/2 
the analog decay b' — ■> t' + W is presumably strongly suppressed by a small mixing 
angle 9' if the mass of the isospin partner t' is larger than my. Once 9' 2 < 1CT 3 
the annihilation decays dominate and single quark decays can be ignored for most 
practical purposes. Similar considerations apply to decays of isosinglet quarks D 
into b + Z,t + W or c + W which are also inhibited by small mixing angles. For the 
subsequent discussion these single quark modes will be ignored. The dominant 
decays are thus all proportional to the wave function squared, the branching 
ratios r a b/^r a fe are independent of the wave function, and as a result can be 
predicted unambiguously. 

The WW mode needs further specification: The decay of rj b > into WW proceeds 
through the virtual isospin partner t' and hence depends on the mass of t'. Present 
limits on the deviation of the p— parameter from one as deduced from electroweak 
precision measurements limit the mass splitting Am between mt> and my in a 
drastic way and Am < 200 GeV can be considered a generous upper bound. For 
simplicity we shall assume equal b' and t' masses ( P qi = 1/4 in eqs. (|5|,|T2"|)) in the 
subsequent numerical evaluation. 
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For bound states t)d from isosinglet quarks the WW decay may in principle 
proceed through virtual quarks of charge 2/3. However, the D-t-W -coupling is 
suppressed by small mixing angles and the WW mode can therefore safely be 
neglected. The HZ mode is absent as a consequence of an = 0. 

QCD corrections to annihilation decays are at present available for hadronic 
(gg—), for 77— and for the ti— decay modes. 
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1 + 4 ^(/ i ) 1 
v 3 7T /3 t 



A(A) + P(A)ln 



l-A 



+ Q(Pt 



where nf — 6 denotes the number of light flavours (including top). In the last 
equation of (|i~3"D the functions A(p t ), P(Pt) and Q(/3 t ) are given by 
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where fit — y 1 — Am 2 /M 2 . The first factor in the corrections for it has its origin 
in the rjQ vertex correction, the second factor is given in ||. Finally, the scale in 
a s is adopted to M in our numerical calculation. The QCD corrections to the 77 
decay rate are negative at the level of about -10% whereas the corrections to the 
hadronic decay increases the decay rate by about +35%. The QCD corrections 
to T t t depend on the mass of the boundstate and vary from +40% for M v = 300 



GeV to -20% for M v = 2000 GeV. 

Since 77, jZ, ZZ and WW can be identified in the limit Mz, Mw « Imyp 
after appropriate adjustment of the gauge coupling, the QCD corrections can be 
assumed to cancel to a large extent in the ratios listed in eqs. |S], |9|). Corrections 
for the HZ mode are missing and we shall simply adopt eq. (0) for the present 
purpose. The normalized decay rates R a b are shown in figs. 4a and b for r\y and 
7]b respectively, the resulting branching ratios V abj^^Y ^ are displayed in figs. 
4c and d. The pattern is significantly different for the two cases: 
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The branching ratio of rjy into HZ is sizable and dominant for a large mass 
range. Also the WW mode is fairly important, amounting to about 1% of the 
hadronic mode, which is presumably difficult to detect experimentally. The 77, ZZ 
and 7Z modes are all relatively small with their ratio basically governed by their 
relative couplings and the statistical factor 1/2 for 77 and ZZ . 

T W w ■ ? zz : IV : T 77 « 1 : 0.462 : 0.064 : 0.021 (17) 

This expectation for asymptotic quark mass values is well met above > 
150 GeV. One thus anticipates a large production rate for HZ and even WW will 
become fairly important. 

This is in contrast to the situation for r] D decays: HZ and WW are absent, 
and the asymptotic relation reads 

T 77 : T yZ : Y zz « 1 : 0.60 : 0.091 (18) 

In this case 77, 7Z and ZZ constitute the dominant visible modes with a branch- 
ing ratio of about 10~ 4 . 

The resulting product of production cross sections times branching ratio for 
SSC and LHC energies are displayed in figs. 5a and b for several interesting r]y 
channels, indicating large rates in particular for HZ. The production and decay 
of t]d is shown in figs. 5c and d indicating that rjo could be observed in the 77, 7Z 
and ZZ modes with characteristic signals. 

To summarize: 

Production and decay rates for bound states of exotic quarks can be predicted 
unambigously. The uncertainty from the QCD potential, a s and the structure 
functions appears to be well under control. The rates for pp — > r]v(-^> HZ) + X 
are particularly promising as a consequence of the large branching ratio. The 
branching rations into WW, ZZ, 7Z, or 77 are significantly smaller but in any 
case larger than the corresponding branching ratios of a toponium resonance. 

Throughout this paper the SQD mode has been ignored completely, a clearly 
oversimplifying assumption. The mixing of both b' and D with other quarks is 
presumably small and it seems difficult to make any hypothesis. Instead we com- 
pare the total annihilation rate (for Vj with = 300 MeV) with the single 
quark decay rate (through the charged current) in fig. 6 for an assumed mixing 
angle of 10~ 2 . This illustrates that annihilation decays would be clearly dominant 
for b' under this assumption. For a complete treatment of rjo we would have to 
include FCNC single quark decays D — > qZ. In view of the completely unknown 
mixing angle we refrain from a detailled analysis. It is evident from fig. 6 that for 
mixing angles of O(10~ 2 ) or less and mo below ~ 500 GeV the main conclusions 
remain unaffected. 
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Figure captions 

Fig. 1 Predictions for R(0) 2 /M 3 of the S-wave ground state function of M 
for the potential Vj with A^L=200, 300 and 500 MeV (dotted, solid and 
dashed lines). 

Fig. 2 Ratio between the radiatively corrected production cross section for 

pp — > rfy d + X and the lowest order result for \/S=16 and 40 TeV. We 

' (4) 

use MT set Bl parton distributions with A-=-^=194 MeV and work in the 

1 MS 

DIS factorization scheme. 

Fig. 3 a) Production cross section for pp — > r\yp + X in NLO for v^S=16 and 40 
TeV. We use MT set Bl parton distributions with A^-=194 MeV and 

1 MS 

work in the DIS factorization scheme (potentials as in fig. 1). 
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b) Dependence of the cross section on the parton distribution functions: 
MT set Bl H (solid), MT set SN |] (dotted), MRS set B200 gO 



(short-dash-dotted), KMRS set B Q (long-dash-dotted) and GRV p 
(dashed). 

Fig. 4 Ratios R a b = r ab /T 77 (ab = HZ,ti,gg, WW, ZZ,jZ) of r)y (fig. a) and tjd 
(fig. b) as functions of M and corresponding branchings ratios (figs, c and 
d) for a top mass of 140 GeV in T^. 

Fig 5. Cross section for r\y (fig. a,b) and rjo (fig- c,d) production including NLO 
corrections multiplied by the branching ratios at = 40 TeV and 16 



MS 

work in the DIS factorization scheme (potentials as in fig. 1) 



TeV. We use MT set Bl parton distributions with A^L=194 MeV and 

^ MS 



Fig 6. Total widths for r\yp from annihilation (see fig. 4) compared with the 
widths from single quark decays into b + W. 
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